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Abstract 

We reformulate time evolution of systems in mixed states in terms of the classical 
observables of correlators using the Weyl correspondence rule. The resulting equation 
of motion for the Wigner functional of the density matrix is found to be of the 
Liouville type. To illustrate the methods developed, we explicitly consider a scalar 
theory with quartic self-interaction and derive the short time behaviour with the 
non-interacting thermal density matrix as initial condition. In the scalar case, the 
complete correlator hierarchy is studied and restrictions are derived for spatially 
homogeneous initial conditions and systems with unbroken symmetry. 
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1 Introduction 



The study of statistical ensembles of quantized fields out of thermal equilibrium has 
a number of interesting applications. Prominent examples are the study of hot QCD 
matter in heavy ion collisions or the behaviour of matter in the course of cosmological 
phase transitions. In contrast to equilibrium theory, where considerable success has 
been achieved in recent years Q, the theoretical framework for the description of 
processes far from equilibrium still lacks a systematic basis. Although some attempts 
have been made to derive the first few stages of a cut BBGKY-hierarchy on the 
basis of equal time correlators H and using other methods ||, a comprehensive 
understanding of the interrelation between observables and the abstract Schrodinger 
evolution for systems of quantized fields is still incomplete. 

Physics of non-equilibrium processes typically involves tree elements. Firstly, the 
description of an initial state by means of a number of - in principle measurable 
- observables. Secondly, time evolution defined by a dynamical equation together 
with the corresponding Hamiltonian operator, and third, the observation of some 
quantities at later times. The second and third step can be discussed on the basis of 
fundamental principles. In particular, time evolution of quantum mechanical systems 
is encoded in the von Neumann equation which determines the time evolution of the 
density matrix operator. Integrating time evolution to later times, it is straight- 
forward to extract expectation values for operators. Only the first step requires 
a statistical assumption about the system, since complete knowledge of the initial 
density matrix is usually not accessible. 

In this paper we focus on the dynamical aspects of statistical ensembles of quan- 
tized fields but reformulate time evolution in terms of classical quantities instead 
of operator evolution. This casts the theory in a form which appears more natural 
when one is interested in observables. Using the effective action as dynamical object 
generating equal time correlators, this problem has been discussed recently Qj. 

Based on the Weyl correspondence principle, we consider the classical counter- 
parts for quantum mechanical operators and the density matrix. Starting with the 
von Neumann equation, we focus on the construction of a dynamical equation of the 
Liouville type for the classical counterpart of the density matrix operator. In order 
to illustrate the abstract result, we derive an expression for the Weyl transform of 
the non-interacting equilibrium density matrix and use it as an initial condition for 
the time evolution in an interacting system. We find that mass resummation is nec- 
essary in order to obtain a sensible physical result. Furthermore, we explicitly derive 
the hierarchy of equations of motion for general initial conditions for the model of 
a scalar field with quartic interaction. Particular conditions on field correlators for 
systems which exhibit spatial homogeneity and with unbroken symmetry are studied 
systematically. 

2 Quantum-mechanics 
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2.1 General 



In order to outline the concepts applied, we demonstrate the basic steps using a zero- 
dimensional quantum mechanical system as example. A harmonic one-dimensional 
oscillator is defined by the Hamiltonian H = ua^a which Fock-space is spanned by 
the eigenstates of H. Systems in mixed states are described by a hermitian density 
matrix operator p, normalized to unity. Expectation values of operators A are given 
by the trace < A >= trpA which is independent of the basis chosen. It can be 
calculated in energy-eigenstates or in position eigenstates as well, 

< A >= ^2 < n\pA\n >= dx < x\pA\x > . (1) 
n J 

2.2 Weyl-correspondence principle 

Each operator is sub-ordinated a classical observable in a state by taking the trace 
of the product with the density operator. An observable defined in such a way 
depends on the particular mixed state of the system under consideration. One may 
thus ask if it is possible to define the corresponding classical quantity independent 
of any particular density matrix, which in a more universal way permits to calculate 
observables in a given state. This is in fact possible by the Weyl correspondence rule 
0. Suppose we can decompose an operator as 



A = J dadTA(a,T)e iap+iT£ , 



where x and p are the canonical position and momentum operators resp., the corre- 
sponding classical function is defined to have the same Fourier transform, 

A{p,x) = J dadTA{a,T)e iap+iTX . (3) 

So far, this is only a definition of the classical counterpart. In order to justify it, we 
consider the quantum-mechanical expectation value of an operator in a mixed state 
and represent it in terms of the classical functions of the operator and the density 
matrix, 

< A >= f dx j dadTda'dr'A(a,r)p(a',T') < x\e iaf,+iT£ e i<7 ' p+iT ' £ \x > . (4) 

Using \p,x] = —i, x\x >= x\x >, exp(irp)|x >= \x — r >, and the Campbell-Baker- 
Hausdorff formula, 

e (A+B) = e A e B e -±{A,B]^ if ^ = ^ = 

one finds 

J dx J dadrda'dr'pia, r)A(a', t') < x + a + a'\x > e \{^-a'r')+i{r+r')x ^ (5) 
and finally on account of the orthogonality of position eigenfunctions, 

< A >= 2vr J dadrp(a, t)I(-ct, -t). (6) 
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That reads, in terms of the classical counterparts of A and the density matrix, 

< A >= J dp J dxp(p,x)A(p,x). (7) 

We conclude: A quantum-mechanical expectation value can be written in a form 
analogous to classical mechanics where the classical counterparts of the operator and 
the density matrix are given by the Weyl correspondence principle. The classical 
density operator p(p, x) is commonly called Wigner function. 

It should be pointed out that the correspondence rule implies a symmetric oper- 
ator ordering prescription. Thus, information about ordering of operators is lost by 
the transformation. However, this is an intrinsic property of any correspondence rule 
since the classical counterparts always commute. In other words, generally, there 
exist distinct operators which correspond to the same classical counterpart. 



2.3 Wigner functions 

Given a density matrix, one may ask how to calculate the classical counterpart. For 
that purpose, we express the Wigner function in terms of the matrix element, 

< z \ p \y >= J dadTp(a,T) < z\e iap+iT % >, (8) 

which evaluates, using exactly the same operator relations as above, to 

<z \p\y>= J drp(y - z,T)e% T{y+z) . (9) 

If we let z = x — 20", y = x + \cr } and perform a Fourier transformation, then 

J da^P <x- ^a\p\x + hj >= J dadTe iap+iTX p(a, r) = p(p, x), (10) 



which is commonly used as definition of Wigner functions ||. The offset ±^<t in 
the definition accounts for the off-diagonal matrix elements. They are interference 
terms which - in contrast to the diagonal ones - cannot be given a probabilistic 
interpretation, but nevertheless do contribute to the classical quantities. 



3 Representation in terms of creation and an- 
nihilation operators 

So far, we worked in position (momentum) representation of quantum mechanics. In 
order to find a link to field theory, it turns out to be useful to go over to a Fock-space 
representation. To further develop our theory, we will need some definitions and 
relations. 
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3.1 Definitions and Relations 

We introduce a, a) related to p, x by a = ^lo/2x + ip/\^2io from which we get 
exp(iap + irx) = exp(aa^ — a* a) with a = {it — uja) / 'y/2uj. The displacement oper- 
ator 

D(a) = e aat ~ a * a = e -5l Q l 2 e Qat e- a * a = e sl a lV a * e aat (11) 
has the following properties {jj. 

• Unitarity: 

D\a)D(a) = 1 , b\a) = D(-a). 

• Shift-operator: 

D^ 1 (a)aD(a) = a + a, b' 1 (a)a)t)(a) = a 1 + a*. 

• The states 

\a >= D(a)\0 >= e"2 |Q| V ^=\n > 

n vn\ 



are over-complete. 

• Product: 

• Identity: 

• Trace: 



< p\ a >= e -¥.n 2 +m+^* 



/2 i I 

a a\a >< a = 1 . 
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trA = — I d 2 a < a\A\a > . 

7T J 

• Multiplication property: 

tr(b(a)D(a')) = 7r5 2 (a + a'). 

• Operator Fourier-decomposition: 

On account of this and the pre-last relation, an observable amounts to the convolution 
of the Fourier-components of density matrix and operator, 



tr(pA) = I d 2 a I d 2 a'A(a)p(a')tr(b(a)b(a')) = I d 2 aA(a)p(-a). (12) 
The Fourier components are calculated using 

1(a) = tr (AD(-a)) . (13) 
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3.2 Time evolution 



Time evolution is given by the von Neumann equation, ip = [H,p], i.e. for p(a) = 
tr(pD(-a)), 

f>(a) = itr(p[H,D(-a)}). (14) 



3.2.1 Example: Free motion 

One finds 

[D(— a), H] = uj{aa^ + a*a + aa*)D{a). 

Weightened values of creation and annihilation operators in (|l4| ) can be written as 
partial derivatives, 

tr (paD(-a)) = ^ - \ap{a), tr (pa!D(-aj) = + \a*p(a) (15) 

and the equation of motion turns into a partial differential equation 

( * dp(a) dp{a) \ 

lp = u) \ a ^- a ^r)- (16) 

Its solution is p(ae~ lujt , a* e lujt ) where we imposed the initial condition p\ t=0 = /5(a). 



4 Field theory 

In order to generalize the setting to field theory, a Lorentz invariant formulation 
is called for. In that context we encounter a fundamental problem. A Lorentz 
invariant definition of an inner product is based on a mass-shell condition. A priori, 
the latter - if it exists at all - is not known in an interacting theory. However, 
from a perturbative point of view, it suffices to know the inner invariant product of 
the non-interacting theory. This product does not pre-assume any deeper dynamic 
property of the full theory, but merely is a necessary starting point of a well defined 
perturbative expansion, which simply absorbs free motion into the interaction picture 
representation. 



4.1 Lorentz invariant product in the free case 

Particles of mass m define two hyperboloids ko = ±uJk with energy iv^ = + m . 
Real on-shell functions on space-time can be Fourier decomposed as 

r(t,x) = J (dp+(k)+dp-(k))e-* k *T(k) (17) 

where dp ± (k) = d 4 kS(k 2 — m 2 )#(±/co) are the measure components invariant un- 
der orthochronical homogeneous Lorentz transformations G C + . An invariant inner 
product for real fields is given by 

(T,£)= f{dp + (k)+dp-{k))T*(k)a(k)= f dp + (k)(r* (k)a(k) + r(k)a* (k)) (18) 
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which reads in position space 

(T,S)=i J d 3 x(T*(t,x)5iS(t,x)-(a t T*(t,x))S(t,x)). (19) 

The direction dt characterizes the space-like hyper-surface on which the quantisation 
is performed. The inner product is constant along that direction. 
The field operator in interaction picture representation 

*/(*,*) = (^372 / d» + (k)(b(k)e- ik * + fct(k) e *») (20) 

decomposed into creation and annihilation operators with commutator 

[6(k),6t(k)] =2w fc <5 3 (k-k / ) (21) 

evolves with the free Hamiltonian 

:H : = J dp+(k)u; k tf(k)b(k). (22) 

The rescaled definition of creation and annihilation operators o(k) = y/2uj k b{k) used 
in the previous section in the zero-dimensional case, and commonly advocated to 
in textbooks, was introduced and discussed by Newton and Wigner ||. One may 
regard the corresponding wave function as a Lorentz-frame dependent probability 
amplitude since their scalar product in position space takes the same form as in the 
Schrodinger theory. 

The exponent of the shift operator evaluates to 

i(E j nj)+i(T,*/) = J dp + (k)(f3{k)b\k)-(3*(k)b(k)), /J(k) = tr(k) - u k a{k), 

(23) 

where we took into account that 11/ = $7. For two reasons, the definition of (3 is 
chosen to differ by a factor \/2uj k from the analog relation for a. The exponent 
formally has a Lorentz invariant form and the properties and relations derived in the 
previous section can be carried over to field theory straightforwardly. In particular, 
the coherent states chosen as basis are created by acting with the shift operator 
functional in interaction picture representation 

Dj[P] = exp || d M +(fc)G9(k)&t(k) - /3*(k)6(k))| (24) 

on the vacuum. Note that only the introduction of a Lorentz invariant product 
renders the exponent time independent, which will turn out crucial in what follows. 

4.2 Full dynamics 

It is useful to split off free (unperturbed) motion by going to the interaction picture 
representation of p which obeys the evolution equation 

ik{t) = [W I {t),p I {t)\. (25) 
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We introduce a Wigner functional decomposition 

pj(t) = Jv 2 p Pl [P](t)DM (26) 
and its inverse, the Wigner transform of the interacting density matrix, 

p I \fl(t)=to(pi(t)D I [-0\), (27) 
which is subject to the equation of motion 

p I W(t) = itr(p I (t)[Wi{t),D I [-0\]). (28) 
The evaluation of time dependent weightened operator expectation values, 

tr (p{t)A) = tr (Mt)ij(t)) = / V 2 (3 Pl [-(3}(t)MP] (29) 
also calls for calculating the interacting operator Fourier components 

A I \p){t)=tv(A I {t)D I [-0\). (30) 

4.3 Equation of motion 

Let us examine the commutator [Wf(t), D{— (3]] at the r.h.s. of the equation of motion 
more closely. We suppose that the interaction is a polynomial in the fields and 
momenta, to wit, also a polynomial in creation and annihilation parts. Since D acts 
as a shift operator, (We drop the index I in what follows.) 

D[-0\b(k) = (6(k) + P(k))D[-0\, D[-/3]&t(k) = (6 f (k) + (3*(k))D[-/3] (31) 

we get for the commutator 

D[-0\] = (W[$] - W[<5> + V])D[-0\ (32) 

with 

*(*,x) = J dp + (k)(f3(k)e-* kx + fT (k)e^). (33) 



In order to draw the operator- valued commutator out of the trace in Eq. (28), it is 
useful to represent b and b^ by suitably chosen functional derivatives acting on D, 

b(k)D[-0\ = (2^^y " \m) D[-fJ\, 
b\k)D[-P] = {-^kj^ ~ i/T(k)) D[-(3}. (34) 

Products of operators b and can be represented by replacing them with their partial 
derivative counterparts in inverse order. The terms with the factor | in front originate 
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from the fundamental commutator and are of relative order fi. The 'classical' field 
operator 

*<*' X > - <2^mJ d " + ^ {e-^-^jki) <35) 

commutes with the displacement ^(t, x) which permits us to write the commutator 
as 

(W[$-ht,]- W [§ + ly]) D [-0\. (36) 



The equation of motion (^8() for the interacting part of the density matrix thus 
assumes the form of a Liouville equation 

= iL[m)p[mi = <w$ - - w$ + \m m 

A formal solution can be constructed by standard methods, 

p[f3](t) =Texpl [ iJ*dt'L[P](t')}plf3}(to) (38) 

using the time ordering operator T . This equation contains the full dynamical infor- 
mation on the Wigner transform of the density operator. The equations of motion of 
the c-number counterpart of p can be obtained by a functional integration with the 
classical field and momentum in the displacement operator similar to Eq. (|3]). 

4.4 Equilibrium distribution 

In order to illustrate the formalism derived, we consider a non-interacting system at 
thermal equilibrium in which interaction is switched on at some initial time. In gen- 
eral, canonical ensembles in thermal equilibrium with temperature T are character- 
ized by the density matrix p eq = e~ H l T . Due to standard arguments in perturbation 
theory, this density matrix can be rewritten as 



p eq = e- H °/ T Ui(to-i/T,t ), U I (t -i/T,t )=Texp^-iJ t ° t dt'Wi(t')j . 

(39) 

The corresponding Wigner function in interaction picture is given by 

Peq [(3](t) = tr (er Ha / T Ui(t - i/T,to)D[-0\) . (40) 

Let us calculate the non-interacting equilibrium Wigner function defined by 

p%[(3]=tr(e- H ^ T D[-(3]) (41) 

which we will assume to be the initial preparation of the system. For that purpose, 
we consider the non-interacting part with one creation operator inserted between 
the statistical factor and the shift operator. Commutating to the right, using the 
cyclicity of the trace and the commutator with Hq we find 

tr (e-^ o/T 6 t (k) J D[-/3]) = -(3*{k)p° eg [(3} + e w(k)/T tr (e^ ^ (k)£> [-/?]) . (42) 
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On the other hand one gets in the trace by acting with the differential operator of 
Eq. d||) on D, 

(~ 2u)k im ~ ^* (k) ) p ° egm = tr ( e ~ Ho/Tbi w D i-®) • ( 43 ) 

Combining this property with the previous equation and functionally integrating with 
respect to /3, one finds, up to a normalisation constant, 

p%\0[ = exp || dp+(k) Q - n-(k)) /3(k)/3*(k)} (44) 

with n ± (k) = (1 - exp(±a; fc /r))- 1 . 



4.5 Non- interacting partition function as initial condi- 
tion 

We consider the non-interacting partition function (O) as initial condition in the 
time evolution (38) and expand the exponent in powers of the coupling constant. To 
be specific, we consider a quartic interaction term in the Hamiltonian, 



W[$\ = ^ J d 3 x^(x). 



(45) 



We act with the differences of Eq. d37l) on p° [(5] and obtain 



{Hx)+ l -^{x))pl q [ 



T(x)p%m, (Hx)-~v(x))p° eq i 



-r*(x)p° eq [P] (46) 



with 



\x) = J dp + {k) (/3(k)n+(k) ( 



+ /?*(k)n"(k)e i 



Green functions arise in the commutator of the differences with T, 



(47) 



G(x-y), 



*(x)±-*(x),r*(y) 



-G*(x-y), (48) 



with 



G(x - y) 



1 



(2tt) 3 

The first order term in A is found to read 

A 



dp + (k)(e- ik< - x -^n-(k) - e ik< - x -^n + (k)). 



L[(3}(xo)p° e M = 77 / d 3 xW[(3}(x)p° e M 



4! 



(49) 



(50) 



with 



W(x) = 6G(0)(r*(x) 2 - T(x) 2 ) + (T*(x) 4 - T(x) 4 ). 



(51) 



For higher orders, it turns out useful to introduce a graphical representation. We 
denote a factor T, T* by an external leg. Green functions and their complex conjugate 
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Figure 1: First order contribution. 





correspond to internal lines with an without star. Each power of the interaction term 
gives rise to a vertex. 

In that way, the second order contribution 

L[(3}(x )L[(3](y )p%i(3] = ^- J d 3 x J d 3 yW[(3](x, y)p%[/3] (52) 

can be represented by the graphs (|2|) trough (|5|), where legs on the l.h.s. and r.h.s. 
correspond to the arguments x and y of T. W[f3](x,y) involves contributions with 
two, four, six and eight external legs and is real. 

The time dependent density matrix in /^-presentation with the given initial con- 
dition is thus given by (We take to = 0.) 

ix /•* 



P [f3](t) = (i + ^J dx J d 3 xW(x)[P]- 



X ~ ft dx f X ° dy I d A x I d d yW(x, y)[(3] I p%\0[. (-v!) 



(4!) 2 7o Vo 



These expressions contain all information about any operator expectation value to 
order A 2 . 
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+ h.c. 



Figure 3: Second order with four external legs. 




+ h.c. 



Figure 4: Second order with six external legs. 
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Figure 5: Second order with eight external legs. 



4.5.1 Particle density 



The extraction of the particle density is achieved as follows. We act with the bilocal 
particle number operator ra(ki,k 2 ) = 6 T (ki)6(k2) in /^-representation on the time 
dependent density matrix fl53|) and set /? = afterwards. In particular, commuting 



n[/?](ka,k 2 ) = 
to the r.h.s. using 

y[0(k),r(x 

6t[/3](k),r*(x 



5f3* (k 2 ) 2 

n+(k) 
"(2vr) 3 / 2 ' 

w~(k) 
~(27r) 3 / 2 ' 



/3(k 2 ) -2 Wfcl 



1 



^(ki 



/3*(ki) 



(54) 



-ikx 



[&[/?] (k),r(s) 



[6[/?](k),r(x)] 



n ( k ) ifc* 

I c 1 



(2vr) 3 /2 

n+(k) 
= (2tt) 3 / 2 ' 



ikx 



(55) 



essentially replaces two external legs by plane waves weightened with distribution 
functions n ± (k). Thus, after having set ,9 = 0, the set of graphs that can contribute 
are those with two external legs only. Moreover, explicit calculation shows that the 
diagrams in Fig. (1) vanish and the first two graphs of Fig. (§) cancel after integration 
over x. The zeroth order term gives (u>i = <^ki) 



25 3 (ki-k 2 Vin+(ki) 



(56) 



and the third and fourth diagram in Fig. (||) contribute, after integration over x, y 
and x ,y 

, G 2 (0) 



A^ d (k! 



(2wi 



r (n + (ki)-rc-(ki))sin^(u;if). 



(57) 



The sunset-diagrams in Fig. (0) are the ones far most complicated. They contribute 
(r = x - yo) 



A 



2 ^ 3 (ki-k 2 ) 
6(2tt) 9 



dxr 



1(x,t) 



■t'o 



rfyo(e~ iLJlT n-(k!) +e^ ir n + (k!)) / d 3 xe iklX / 3 (x, r) + h.c, 



^(e~^ T+ipx n-(p) - e^ r ~ ipx n + (p)). 

2Wr, 



(58) 



Let us discuss the contribution (|57| ) in more detail. It exhibits a harmonic time de- 
pendence and does not decrease with time as one may expect for the particle number 
density which presumably approaches an equilibrium value. The situation is even 
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Figure 6: Resummation of mass tadpoles 



worse for off-diagonal operator expectation-values. One finds for the contribution 
from the third and fourth diagram in Fig. (|2|) to the operator 6(ki)6(k2) 

X 2 5 3 (ki + k 2 ) G ^ fsin 2 ^it + i(tu\ - sinc^t coswii)) (n + (ki) - rc~(ki)), (59) 
(2u;i)' ;i v / 

which even contains terms linear in time. This indicates a breakdown of the expansion 
of the time ordered exponent in the evolution (|38|), A crude estimate of the domain 
of validity of the expansion performed can be given by the requirement that the 
exponent remains smaller than unity, i.e. for times Xt <C 1. 



4.6 Resummation 

As well known in equilibrium finite temperature field theory, consistency in the in- 
frared regime requires a resummation to be carried out to reorganize the pertur- 
bative series |IJ|. As anticipated in the example given, for non-equilibrium density 
distributions the necessity of resummation also arises for dynamical reasons. Here, 
resummation boils down to add a mass square M 2 to the original Hamiltonian and 
compensate it by a two-vertex 



-J d 3 xM 2 $ 2 (x) (60) 



in the interaction part. This term changes the Liouville-operator in the equations of 
motion, which gives rise to a number of additional diagrams. Those are such that 
each graph with closed loops oc G(0) gets a mass counterpart in a way that they 
compensate for the choice 

M 2 = ^G(O), (61) 

as illustrated in Fig. (||). 

The resummed mass m 2 + M 2 appears now hidden in the Lorentz-covariant vol- 
ume element d/j, + and the corresponding energy eigenvalues \/m 2 + M 2 + k 2 both 
appearing in G(0). In that sense, Eq. ( |6l|) is a consistency condition for the mass 
counter term. It is interesting to observe that not the first order dynamical con- 
tributions to the particle number, which vanish, enforce the resummation, but the 
time evolution of the second order terms ([57]) necessitate the reorganisation of the 
perturbative series. The redefined theory gets its second order contributions (|58|) 
solely from the non-local sunset-graphs. 
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5 General initial conditions 



We generalize the particular density matrix considered in the previous section by a 
Wigner density functional of the connected graphs expressed in a functional power 
series in f3 (See Appendix.), i.e. 



logp[/3](i) 



oo 

EE 

n=l {a} 



i (n / d^ik^j /W(t; k ls . . . , k„)/?^ (kx) . . . /T"(k n ), 



(62) 

where elements of the set {a} = {a±, . . . , a n } can assume values + and — which refer 
to starred an unstarred /3's respectively. The correlators f n may also explicitly depend 
on time. For a given order in f3, not all functions f n are independent. They have to 
be symmetric with respect to pairwise exchange of two momenta, both corresponding 
to starred {a = +) or unstarred {a = — ) pairs of /3's. Furthermore, since physical 
density matrices have to be hermitian, the relation D[{3] = D[— j3y implies invariance 
of p\(3\ with respect to simultaneous hermitian conjugation and sign flip of (3. This 
implies for the correlators / 



ft } \t-k 1 ,...,k n ) = (-) n ft } (t-M,...,k n 



(63) 



The number of independent correlators is (n + 2)/2 for n even and (n + l)/2 for n 
odd. 

Let us study the r.h.s. of the equation of motion (|37|). Acting with the Liouville 
operator on p{(3] = exp(log p[(3]) gives 



L[(3]p[(3] 



d 3 x^(x) (3[4>(x), [$(x),log / o]][4>(x),logp] + 

1 



+([<$>(x), log p\y + [*(ar), [*(a;), [$(x),logp]]] + -^(x)[$(x), log p]\ p. (64) 

The first few terms may best be studied by introducing a graphical notation. We 
represent log p by the graphs depicted in Fig. ([?]) . 

The blob in the middle denotes the correlators /, and each external leg a factor 
containing (3 or (3* . We do not explicitly split the sum over {a} into their summands 
- each graph is understood as the sum over all independent f n 's for a fixed order 
n. Acting with <&(x) on an external leg cuts it off and replaces it by a plane wave. 
We denote cut legs by a cross. The factor ^ is represented by a cut leg without 
blob since it corresponds to a unit cut one-point correlator. The integration over 
space graphically joins all crossed legs in the point x. The time derivative idt in the 
equations of motion is symbolized by a dot. The motion of the system to order /3 4 is 
depicted in Figs. (|8|-|Tl|)- 



Let us clarify the magnitude of these graphs in orders of %. The relevant quantities 
have the following dimension. 



object 


L[(3] 






P 


5/5(3 


fn 


order in h 


1 


1/2 


3/2 


-1/2 


-5/2 


-3n/2 
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Figure 7: n-th order term in log p. 




Figure 8: First stage of the hierarchy of the equations of motion. 




Figure 9: Second stage of the hierarchy of the equations of motion. 
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Figure 11: Fourth stage of the hierarchy of the equations of motion. 
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Figure 12: n-th order stage of the equations of motion. 



One can easily check that all graphs are of the same order in fi except for the last 
in Fig. (|To|) and the third at the r.h.s. in Fig. (11), which are suppressed by a factor 



Ti . In general, all those terms originate from the last summand oc Vf 2 in Eq. (|64|). 
A systematic expansion thus requires to keep all but those terms of higher order. 
However, the structure of the the whole hierarchy is such that the time derivative of 
correlators of order n involves correlators of order n + 2 and lower on the r.h.s. in the 
equations of motion. The Ti 2 suppressed terms are of order n — 2 in 0. In Fig. ([1^), 
for any n, all graphs are shown which involve higher or equal order correlators at the 
r.h.s. They are responsible for the fact that cutting the hierarchy gives a non-closed 
system of equations. Even worse, the hierarchy cannot be cut at some order by the 
condition f n = for n grather than some fixed number, since that would contradict 
the equations for /„ containing nonvanishing contributions at the r.h.s coming from 
lower order correlators. 



6 Particular configurations 

Many systems of physical interest exhibit particular symmetries which imply restric- 
tive conditions on the correlators f n . 
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6.1 Spatially homogeneous systems 



Spatially homogeneous systems are characterized by the invariance of operator prod- 
ucts with respect to translation in space, to wit, 

tr ($(t, xi)$(t, x 2 ) . . . x n )p) = tr 0)$(i, x 2 - Xi ) . . . $(x n - Xi)p) . (65) 

The corresponding expression in terms of Wigner transforms is found by replacing 
$(i,x) by <5(i, x) — |\I/(i,x), the density operator by its transform and the trace by 
a functional integral. Let us pick that term which only contains Xj)'s, act with it 
on p[(5] and set (3 = afterwards. Without loss of generality we may furthermore set 
t = since time evolution commutes with spatial translations. The condition (|6^) 
implies 

(ft / dn+ih)^^ . . . ,k n ) = 

(ft / d » + ( k i)^J exp ^fk 1 (x s - xOaJ ftH^i, ■ ■ ■ ,K). (66) 

However, this equation can only hold for correlators which are such that the exponents 
at both sides are equal. That implies the particular form 

/i CT >(k!, . . . , k n ) = /i CT} ' w (k!, . . . , k n )5 (£ vikij (67) 

where the homogeneous correlator includes one redundant momentum. The delta- 
function of the a-averaged momenta generalizes momentum conservation for spa- 
tially homogeneous systems to the correlators. Spatial homogeneity and isotropy is 
a constant of motion, to wit, homogeneous initial conditions evolve to homogeneous 
correlators. 



6.2 Unbroken symmetry. 

Another physically interesting particular situation are configurations with unbroken 
symmetry. They are characterized by vanishing expectation values of the field oper- 
ator, i.e. 

tr ($(t,x)p\ = 0. (68) 

Again, this condition implies a restriction on the Wigner transforms. The operator 
$(i,x) — ^^(t,x) decreases (<J>) resp. increases the number of /?'s in the trace by 
one. Comparing all terms with equal number of /3's, which have to vanish separately, 
one easily finds that all odd correlators have to vanish, 

fW> hom (k 1 ,...,k n ) = 0, nelN oM . (69) 

Closer inspection of the graphical representations given in Figs. (|8]-[lT|) and their 
higher order counterparts reveals that the time derivative of odd correlators also 
vanishes if they represent a system with unbroken symmetry, which is necessary for 
consistency. The symmetry conservation is a consequence of the invariance of the 
Hamiltonian under sign flip of the field. 
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7 Conclusion and Outlook 



In this paper we reformulated the von Neumann equation for the density matrix in 
terms of the classical observables of field correlators using the Weyl correspondence 
principle and the language of Wigner functions. We found a Liouville equation for 
the Weyl-transform which expands into an infinite hierarchy of equations of which we 
showed that they cannot be cut consistently at a finite order. Time evolution is known 
exactly, and can be - at least perturbatively - be integrated to later times. Further 
modelling of physical scenarios, however, includes additional statistical assumptions 
entering in realistic initial conditions on the density matrix. Based on the present 
work, investigations in that direction are currently under progress. 
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Appendix: Proper and Improper Correlators 



We call the set of correlators / n (t;ki, . . . ,k n ) proper correlators if they cannot be 
written as a product of two functions, i.e. 

fn ^(^! ki, ■ ■ ■ j k n ) ^ u m (t] k ri , . . . , ls. rm )v n — m {t\ k rm+1 , . . . k rn ) (70) 

where m < n and {r%, . . . ,r n } denotes an arbitrary permutation of the set of indices 
{1, ...,n}. In particular, proper correlators are not the product of correlators of 
lower order. Conversely, improper correlators can always be written as a product. 

We now prove that the exponentiated logarithm in the definition Eq.(|62|) exactly 
produces all improper correlators which are necessary to describe a general density 
matrix p[0\. Let us for the moment drop the plus/minus index referring to starred 
and unstarred components and consider all contributions to a given order n. (Here 
order is meant to be the number of (3 and (3*.) The possible types of improper 
graphs may be described by a partition of n in the following way. Be {h, ■ ■ ■ ,l n } 
a set of natural numbers l u G IN o such that J2u=i lv v = n i then this set describes 
an ensemble of proper correlators which contains l v times - l v may also be zero - a 
proper correlator of order v, and which altogether represents one improper correlator 
of total order n. We determine the combinatorial factor of that improper correlator. 
It has n momenta which can be attached in n! permutations to the /3's. However, 
in this way we over-counted indistinguishable contractions among the momenta of 
the proper subgraphs. For each subgraph, we have to divide by the number v\ of 
permutations within a subgraph. If {1} contains a subgraph of order v more than 
once, l v > 1, permutations of those are over-counted by a factor l v \. Finally, we sum 
over all possible partitions {/}, which gives the n-th order summand of p[(3\ 

^ r 1 / 



M {1} 



One can convince himself that restoring the ± index contained in the sum over {a} 
does not change the symmetry factors. The complete (3 transform is given by the 
sum 

OO -i 

where po[P] = 1 encodes unit normalisation of the density operator. The factor 1/n! 
has to be inserted to be consistent with the definition (^) which also involves this 
factor for each correlator. 



In order to show that the exponentiated expression (62) equals that expression, 
we write the exponentiated sum in the r.h.s as a product of single exponentials. 
Then, the i-th summand in the expanded exponentiated n-th term in the r.h.s of 
( p^ ) exactly corresponds to the term with /„ = i. The correct factor (n!)*i! in the 
denominator comes from the i-th power in the Taylor-expansion of the exponential 
function. The sum over ,l n } just amounts to collect all terms in the product 

of the expanded exponentials which belong to order n. This completes the proof. 
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